Optical Line Emission in Brightest Cluster Galaxies at < z < 
0.6: Evidence for a Lack of Strong Cool Cores 3.5 Gyr Ago? 

Michael McDonald 1 - 2 
ABSTRACT 

In recent years the number of known galaxy clusters beyond z > 0.2 has 
increased drastically, with the release of multiple catalogs containing >30,000 
optically-detected galaxy clusters over the range < z < 0.6. Combining these 
catalogs with the availability of optical spectroscopy of the brightest cluster 
galaxy from the Sloan Digital Sky Survey allows for the evolution of optical 
emission-line nebulae in cluster cores to be quantified. For the first time, the 
continuous evolution of optical line emission in brightest cluster galaxies over the 
range < z < 0.6 is determined. A minimum in the fraction of BCGs with optical 
line emission is found at z ~ 0.3, suggesting that complex, filamentary emission 
in systems such as Perseus A are a recent phenomenon. Evidence for an upturn 
in the number of strongly-emitting systems is reported beyond z > 0.3, hinting 
at an earlier epoch of strong cooling. We compare the evolution of emission line 
nebulae to the X-ray-derived cool core fraction from the literature over the same 
redshift range and find overall agreement, with the exception that an upturn in 
the strong CC fraction is not observed at z > 0.3. The overall agreement between 
the evolution of cool cores and optical line emission at low redshift suggests that 
emission-line surveys of galaxy clusters may provide an efficient method of in- 
directly probing the evolution of cool cores and, thus, provide insights into the 
balance of heating and cooling processes at early cosmic times. 

Subject headings: 

1. Introduction 

The high densities and low temperature of the intracluster medium (hereafter ICM) 
in the cores of some galaxy clusters, dubbed "cool cores" (hereafter CCs) due to the steep 
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rise of the density profile and drop in the temperature profile, im ply massive a mounts of 
gas cooling radiatively out of the ICM. Early studies (see review by lFabianlll994l ) suggested 
that cooling flows should deposit on the order of 100-1000 M yr~ x of cold gas onto the 
brightest cluster galaxy (hereafter BCG) in the cluster core. The lack of evidence for such 
vast amounts of cold molecular gas and young stars is often referred to as the "cooling flow 
problem" . In recent years, the understanding of cooling flows has benefited immensely from 
several multi-wavelength observing campaigns, which suggest that ~90% of the cooling is 
halted by some form of feedback, while a smaller but still significant supply of cooling gas 
makes its way to the cluster core. 

The presence of optical emission-line nebulae has b een noted in a large number of CC 



clust e rs to date, the m ost remarkable being Perseus A fjCowie et al.lll980t IConselice et al. 



200 ll ; lHatch et al.l 120061 ). While the ionization mechanism of these filaments remains a topic 



of rich debate, the cool gas i tself most likely origina ted in the hot ICM (e.g., iHeckman et al. 



19891 ; I Crawford et al.l 120051 ; iMcDonald et al.l 120101 ) and is either undergoing initial inf all or 
has been uplif ted from the core ei ther by a buoyant bubble (e.g.. IChurazov et al.ll200lf ) or a 
radio jet (e.g.. IWerner et al.ll201ll ). An important aspect of these nebulae which has received 
little attention is their evolution over time. The fraction of galaxy clusters with optically- 
emitting BCGs as a function of redshift may provide valuable insights into the relative 
effectiveness of heating and cooling processes in the ICM over cosmic time, which is critical 
for models of galaxy cluster and struct ure formation. . Furthermore, the fact that CC s tend to 



have excess optical line emission (e.g.. lCavagnolo et al.ll2009t IMcDonald et al.ll2010l ) suggests 



that tracing the evolution of optically-emitting BCGs may indirectly probe the evolution of 
CCs in galaxy clusters, which proves exceedingly difficult using X-ray techniques at high 
redshift. 

The purpose of this Letter is to examine the evolution of optical emission-line nebulae 
in galaxy clusters as a function of redshift. In §2, the data and analysis techniques are 
presented, with a focus on aperture corrections. The evolution of optical line-emitting BCGs 
is presented in §3 over the range < z < 0.6 and the implications of this evolution are 
discussed. In §4 we compare this evolution to that of the CC evolution from the literature 
and speculate on the ability of optical surveys to identify CCs. Finally, in §5 the implications 
of these results and future avenues of research are summarized. 



Throughout this Letter, we assume H = 71 km s 1 Mpc 1 , VLm = 0.73, and Q t 



0.27. 



- 3 - 



Sample 


Number of 


Redshift 


Detection 


Reference 


Name 


Galaxies 


Range 


Method 




Wen+09 


39,716 


0.05 < z < 0.6 


• Friend-of-friend algorithm 


Wen et al. (2009) 








• N > 8 




GMBCG 


>55,000 


0.1 < z < 0.55 


• Red sequence matching 


Hao et al. (20101 


Szabo+11 


69,173 


0.045 < z < 0.78 


• Adaptive matcher filter 


Szabo et al. (2011) 



in position, magnitude, 
and redshift 



Table 1: Summary of optically-selected galaxy cluster catalogs used in this study. Further 
explanation of the detection algorithms and sample properties can be found in the individual 
catalogue references. 



2. SDSS Data and Analysis 



In recent years, the number of known galaxy clusters in the range 0.1 < z < 0.6 
has increased significantly due to the application of complex cluster-finding algorithms to 
the wealth of publicly -available data in the Sloan Digital Sky Surve}0 (hereafter SDSS; 



Abazajian et al.ll2009f ). In this study, we consider three recently- released catalogs, each of 
which identify galaxy clusters based on SDSS optical photometry, using different algorithms, 
as described in Table [TJ From each of these catalogs, the BCG positions were cross-referenced 
with the MPA-JHU@ release of SDSS DR7 spectrum measurements, to obtain Ha, H/3, [O II], 
and [O III] emission line luminosities. All spectroscopic measuremen ts from the MPA -JHU 



release have already been corrected for Galactic reddening, following lO'Donnelll (119941 ) . 



In order to compare spectroscopic quantities over a range of redshift, a suitable aperture 
correction must be applied to account for the increased apparent diameter of the 3" SDSS 
fiber aperture at larger distance. In this work, we consider three separate aperture corrections 
in an attempt to determine the total emission-line luminosity from BCGs over the redshift 
range < z < 0.6. These three methods are: 



1. Universal L# Q (r) profile - Using deep, high spatial resolution Ha im aging of four 



BCGs with th e most luminous, extended, optically-emitting filaments from lMcDonald et al. 



( 120101 . l2011al ). we derive a Universal L# Q (r) profile, as shown in Figure [Jji. This is used 
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to estimate the fraction of the emission line flux which falls outside the 3" aperture as 
a function of redshift (or apparent aperture size). 

2. No evolution in < Ljjp > - In Figure [lb the H/3 luminosity is shown as a function of 
angular diameter distance for 12,230 BCGs in the GMBCG catalog. The trend towards 
higher emission-line luminosity at large distance is most likely a result of increasing 
aperture size. Thus, we make the assumption that the mean H/3 luminosity should 
be independent of redshift and simply fit a power law to these data. Applying this 
correction results in a flat distribution of Lh/3 with distance. 

3. LiHa oc L con tinuum - By considering the equivalent width (line luminosity normalized 
to continuum luminosity), we make the assumption that the warm line-emitting gas 
has a similar radial distribution as the stellar content. This appears to be a reasonable 
assumption for an ensemble, as the equivalent width shows no strong evolution with 
angular diameter distance (Figure [Tb). 
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Fig. 1. — Graphical depiction of three separate aperture corrections which have been used in 
this work. Left: Simulated fraction of Ha luminosity within a 3" aperture as a function of an- 
gular diameter distance, Da, for systems with complex, extended Ha filaments: Abell 0478, 
Abell 1795, Abell 2597, Sersic 159-03. The red, dashed line is a fit to these points and rep- 
resents the first of three aperture corrections used in this work. Middle: H^ luminosity as a 
function of Da for 12,230 BCGs in the GMBCG catalog. In order to correct for the increas- 
ing aperture size, we assume that the average Hg luminosity is redshift-independent. The 
dashed red line represents the aperture correction required to remove the observed non-zero 
slope. Right: H^ equivalent width, EWnp as a function of Da- The use of EWup assumes 
that the radial distribution of H/3 luminosity is proportional to the continuum emission. 



These three methods of estimating the total emission-line luminosity based on the central 
3" provide a range of aperture corrections which can then be compared. Both the magnitude 
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Fig. 2. — Total Ha lum inosity within a 3 " diameter ap erture versus redshift for a sample of 
13 galaxy clusters from McDonald et al.l ( 120101 . 1201 lal lbh . Systems with complex, extended 
emission and nuclear emission are depicted with red crosses and blue circles, respectively. 
Horizontal dashed lines represent luminosity limits of 10 40 ' 2 and 10 40,8 erg s _1 , chosen to 
separate extended from compact emission and to allow detection of weak emitters at z = 0.6. 
The curved dotted line represents the estimated sensitivity limits, showing that beyond 
z > 0.4, we may be missing a fraction of the weakly-emitting systems. For z > 0.3, the 3" 
aperture encloses enough of the total Ha luminosity to reliably separate these luminosity 
classes. The upper axis shows the angular diameter of a 3" aperture for a given redshift. 
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and uncertainty of these aperture corrections is largest for nearby BCGs, where the 3" 
aperture encloses only a small fraction of the total light. In Figure [21 the Ha luminosity 
within a simulated 3" apertur e is shown over the range < z < 1 for 13 BCGs from 
McDonald et al.l (120101 . l2011al ). This figure demonstrates that, for z > 0.3, a sufficient 



fraction of the total Ha luminosity is captured by a 3" fiber to differentiate between strong 
(extended, complex filaments) and weak (nuclear, slightly extended) Ha emission. Based on 
Figure |2j we choose limiting Ha luminosities of 10 402 erg s _1 and 10 40 ' 8 erg s _1 to identify 
weak and strong optical emission 

In order to increase signal-to-noise, emission from the Ha, H/3, and [O II] lines were 
combined wherever pos sible. This is don e assuming typical conversions of Ha/H/3 = 2.85 
(case B recombination; lOsterbrockl Il989l ) and [O II]/H/3 = 1.4, which is typical of an H II 
region. For each cluster, given a set of line fluxes and uncertainties, the probability that the 
optical emission lies above the aforementioned luminosity thresholds is computed assuming 
Gaussian statistics. The sum of these probabilities for all clusters within a given redshift bin 
then represents the fraction of clusters with optical line emission, NoLE/N to tai- In order to 
account for contamination due to AGN, we remove all systems with [O III] /H/3 > 3, which 



is atypical for the optical filaments typically found in CC clusters (j Crawford et al.l 11999 
McDonald et al. in prep). 



3. The Evolution of Optical Line-Emitting BCGs from z = to z = 0.6 

In Figure |3j the fraction of galaxy clusters hosting an optical line-emitting BCG as 
a function of redshift is shown. For each of the three catalogs described in Table [TJ we 
show the fraction of clusters with optical line emission in the BCGs, as well as those with 
strong (JjHa > 10 40 8 erg s -1 ) and weak (10 40 2 erg s _1 < Ljj a < 10 40 ' 8 erg s _1 ) emission 
separately. In general, there is little difference between the three cluster catalogs, despite 
the different cluster identification algorithms. All nine panels of Figure [3] show the same 
overall trend, regardless of the emission strength, sample selection, or aperture correction 
used. The fraction of clusters with optical line-emitting BCGs decreases with redshift until 
it reaches a minimum at z ~ 0.3. For z > 0.3, the fraction of clusters with weak line emission 
continues to decrease, while the fraction with strong emission increases. The net effect of 
this opposite evolution is that the total fraction of line-emitting clusters remains constant 
over the range 0.3 < z < 0.6. We caution that, while the line luminosity thresholds have 
been chosen to reduce sensitivity biases, the fraction of clusters with weak optical emission 
may still be under-estimated at high redshift. 

Figure [3] provides strong evidence for a minimum in the evolution of optical nebulae 
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Fig. 3. — Fraction of galaxy clusters with optical line emission (OLE) in the BCG as a 
function of redshift for three different optically-selected samples of galaxy clusters. The 
three samples, described in §2, detect clusters using different techniques and yield > 10, 000 
BCG spectra each. Point type and color represent the three aperture corrections described 
in §2, while solid black points represent data for z > 0.3 clusters with no aperture correction. 
The different catalogs and aperture corrections all yield similar results: there is a decrease 
in the number of galaxy clusters with emission line nebulae from z = to z = 0.3, and an 
upturn in the fraction of systems with strong emission at z > 0.3. 
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in galaxy clusters at z ~ 0.3, or roughly 3.5 Gyr ago. The agreement between different 
indicators (Ha, H/3, [O II]), and different catalogs suggests that selection effects and mea- 
surement errors are minimal. What follows is a brief discussion of the two regimes bounding 
the minimum at z ~ 0.3, and what processes may govern the evolution of optical nebulae 
at these epochs, followed by a summary of the different biases and systematic uncertainties 
which could conspire to produce the observed evolution. 



3.1. Optical Line-Emitting BCGs at z < 0.3 



The evidence for a decline in the fraction of line-emitting BCGs over the past ~3 . 5 Gyr 
is compelling. Interestingly, this timescale agrees with the finding by lMcDonald et al.l (120101 ) 
that Ha filaments extend to radii at which the ICM cooling time is ~ 3-5 Gyr, but never 
beyond. The presence of an app arent equil ibrium between radio feedback from the central 
AGN and the cooling ICM (e.g., ISunl 120091 ) also suggests that cooling has been ongoing for 
a substantial amount of time. Combined, these results suggest that the majority of galaxy 
clusters have been relatively undisturbed over the past few Gyr, during which time the ICM 
has been allowed to cool, resulting in the formation of young stars, a multiphase ICM, and 
the establishment of a feedback cycle between the AGN and the cooling ICM. 

One mechanism which could produce the observed evolution in Figure [3] over the range 
< z < 0.3 is cluster mergers. Based on simulations, iGenel et al.l ( 120091 ) predict that the 
fraction of massiv e halos undergo i ng ma jor mergers drops by a factor of ~4 between z = 1 
and z — 0, while iFakhouri et al.l (120101 ) predict that 50 % of massive halos at z = had 



ZuHond ( 1201 ll ) has shown that a merger between 



their last major merger at z >1. Further, 
two galaxy clusters can fully disrupt a CC, increasing the central cooling time to t coo i ~ 
10 Gyr. This would cut off the supply of cool gas, quenching star formation and resulting 
in a general lack of optical line emission on timescales of ~100 Myr. The ability of major 
mergers to prevent ICM cooling, coupled with the rapid decrease in merger rates for massive 
halos at z < 1 could produce the observed minimum in the fraction of galaxy clusters with 
emission-line nebulae at z ~ 0.3. In the absence of strong mergers at z < 0.3, clusters may 
cool unimpeded until the onset of AGN feedback, producing a significant reservoir of cool 
gas and, thus, a steep rise in the fraction of optically-emitting BCGs. 

Strong AGN feedback is another way of reducing the amount of optical emission at z ~ 
0.3. During a strong outburst, the radio luminosity of the central AGN can surpass the pdV 
work required to remove the CC (e.g., ISunl 120091 ; iMcDonald et al.l 1201 laf) . This imbalance 
between feedback and cooling is seen in a small number of z ~ galaxy clusters, including 
Abell 2052 and Hydra A, which exhibit massive cavities in the ICM. If such outbursts were 



-9 - 



more common in the early Universe (z > 0.3), then ICM cooling would be suppressed. Once 
the CC is destroyed, the AGN would lose its fuel supply and ICM cooling would proceed 
unimpeded, producing the cool gas required to generate significant emission line nebulae at 
z < 0.3. 



3.2. Optical Line-Emitting BCGs at z > 0.3 

In Figure [31 we show that the fraction of BCGs with strong optical line emission increases 
with redshift beyond z = 0.3, while the fraction with weak emission decreases over the 
same range. This decrease in the fraction of weak emitters may be a result of sensitivity 
limitations, despite efforts to set the threshold luminosity above than these limits (Figure [2]). 
The rise in the fraction of strong optical emitters is statistically significant at the 72% and 
>99% confidence levels for the GMBCG and Wen+09/Szabo+ll samples, respectively. This 
upturn suggests a significant amount of cool gas in clusters cores at z ~ 1, roughly 6 Gyr 
after the Big Bang. This second maximum may correspond to the initial inflow of gas from 
the intergalactic medium onto the giant elliptical. Without a dense ICM at high redshift , 



mergers between galaxy groups will be less efficient in halting cooling. Indeed, lEvrardl (119901 ) 
predicts that the majority of the ICM is formed around z ~ 1 in a typical cluster and that 
most of this gas has kT < 1 keV for z > 1. While these results are intriguing, deeper optical 
emission-line surveys are needed at z > 0.3 in order to confirm the observed trends. 



3.3. Biases and Systematic Uncertainties 

The effects of various biases and systematic uncertainties were examined in order to 
assess their impact on the evolution presented in Figure [3j BCG misidentification (~10% 
in CCs; Leisman et al. in prep), bias towards rich clusters at high-z in optically-selected 
catalogs, and decreased sensitivity to Ha at z ~ 0.3 were found to have only a small impact 
on the observed evolutionary trends. As discussed previously, the evolution of optical line- 
emitting BCGs for z < 0.3 is strongly dependent on the choice of aperture correction, but 
the fact that all three aperture corrections agree provides evidence that they are roughly 
correct. For z > 0.3, the effects of aperture correction are negligible, but sensitivity biases 
may lead to an under-estimated fraction of weakly-emitting BCGs. 
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4. Insights into the Evolution of Cool Cores? 



In recent years, the fraction of galaxy clusters hosting CCs has been quanti fi ed vi a 
a variety of methods , as reported in recen t ind ependent surveys by iBauer et al.l ( 120051 ) , 
Vikhlinin et aD fl2007h . febeling et all J2007t ). and ISantos et all tooA I2OI0I) . These studies 
are all based on X-ray data and rely on either the central cooling time or the "cuspiness" 
of the X-ray surface brightness profile to identify CCs. While this is the most fundamental 
way of identifying a CC, which by definition has high density and short cooling times, it is 
also highly inefficient at larger redshift where the number of X-ray counts within the central 
~50 kpc is generally too low to p roperly constrain the ce ntral cooling time. However, a more 
promising method may exist. In iMcDonald et al.l ( 120101 ) we showed that the presence of Ha 
filaments is evidence for cooling flows, with the Ha luminosity being correlated (albeit with 
significant scatter) with both the X-ray cooling rate and central entropy. With this in mind, 
it is relevant to ask the question: Does the evolution of optical emission in galaxy clusters 
trace the evolution of CCs? 

Figure H] shows the evolution of the CC fraction as a function of redshift. The four afore- 
mentioned X-ray studies provide estimates of both the weak and strong CC f raction at various 
epoch s, but the data remain sparse. We turn to the ACCEPT@ database ( jCavagnolo et al.l 



20091 ) to add an additional 241 galaxy clusters with uniformly analyzed X-ray data. Using 



the readily-available central entropy {K ) measurements from the ACCEPT database, the 
CC fraction is estimated over the range < z < 0.4 for weak (20 < Kq < 70 keV cm 2 ), and 
strong (Ko < 20 keV cm 2 ) CCs. This compilation represents the most de tailed picture of 
the evolution of CCs in galaxy clusters to date. In IMcDonald et al.l ( 120101 ) we showed that 
galaxy clusters with extended, filamentary Ha emission also had the highest cooling rate and 
lowest central entropy. Thus, for a large enough sample, the evolution of strong optical line 
emission should mirror the evolution of strong CCs, and likewise for weak CCs. Overlaid on 
the X-ray data in Figure HJ we show the evolution of optical line-emitting BCGs from Figure 
[3j for the GMBCG catalog. We find, overall, qualitative agreement between the evolution 
of optical nebulae and CCs in galaxy clusters, with the exception of any statistically signifi- 
cant increase in strong CCs at large redshift, providing additional evidence for a connection 
between the warm and hot ICM phase and justifying the use of optical line emission as a 
probe of CCs. 

If, in an ensemble of clusters, the presence of optical line emission is indeed a tracer of 
CCs (as Figure H] seems to indicate), the discussion in §3 can apply to CCs as well - namely 
that there is a distinct lack of CCs 3.5 Gyr ago that may be due to increased merger or AGN 



3 http : //www . pa . msu . edu/ astro/MC2/ accept/ 



- 11 - 





u 


Q 

o 


c 






_o 






—i — ' 

u 


U 


D 


o 






1_ 






Li_ 


n 


4 


O 






o 








n 


7 
z 













1 







8 


c 






_o 






-t— * 
o 





6 


o 












L_ 





4 


O 






O 











2 








Santos + 08 
Ebeling + 07 
Vikhlinin + 07 
Bauer+05 




a 

O 
< 



o 
o 

z^ 
□ 

(D 



o 
o 

c 
o 

00 



Fig. 4. — The fraction of galaxy clusters with cool cores (CCs) as a function of redshift, based 
on the literature. The upper panel shows the general evolution of CCs, while the lower panels 
separate into weak and strong CCs. These classifications are based on central cooling time, 
X-ray surface brightness cuspiness, and central entropy. The shaded gray regions represent 
the evolution of optical line-emitting BCGs from Figure |3] (GMBCG catalog). The overall 
agreement between the CC and emission-line nebulae evolution suggests that large surveys 
for optical line emission in cluster cores may be able to accurately probe the CC evolution. 
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activity in the early Universe. 

5. Summary 

The presence of optical emission- line nebulae in the cores of galaxy clusters is often taken 
as evidence for ongoing cooling of the intracluster medium (ICM). Thus, understanding the 
evolution of this warm phase is an important step towards understanding the heating/ cooling 
balance in cluster cores over cosmic time. In this Letter, we exploit the availability of large, 
optically-selected catalogs of galaxy clusters with readily-available optical spectra in the 
Sloan Digital Sky Survey. We find evidence for a decline in the fraction of brightest cluster 
galaxies (BCGs) with optical emission-line nebulae from z = to z = 0.3, suggesting that 
cooling has only recently begun to dominate feedback in cluster cores. For z > 0.3, the 
fraction of clusters with strong optical emission rises, hinting at an earlier epoch of strong 
cooling at z ~ 1. The evolution of optical emission in BCGs is compared to estimates of 
the cool core (CC) fraction and found to agree overall, with the exception that an upturn in 
the strong CC fraction is not observed. This agreement suggests that narrow-band optical 
surveys of high- ,2 galaxy clusters may provide a new, efficient way of probing the evolution of 
CCs at redshifts where X-ray methods become too expensive. In the future, deep, narrow- 
band Ha imaging, which does not require complicated aperture corrections, of both low-z and 
high-z cluster samples will allow for an even more accurate analysis of high-z emission-line 
nebulae. 
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